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Short title: The Complex Geometry of Piecewise Solutions 



Abstract 

An extension of the algebraic-geometric method for nonlinear integrable PDE's is shown to 
lead to new piecewise smooth weak solutions of a class of A^-component systems of nonlinear 
evolution equations. This class includes, among others, equations from the Dym and shallow 
water equation hierarchies. The main goal of the paper is to give explicit theta-functional 
solutions of these nonlinear PDE's, which are associated to nonlinear subvarieties of hypercUiptic 
Jacobians. 

The main results of the present paper are twofold. First, we exhibit some of the special 
features of integrable PDE's that admit piecewise smooth weak solutions, which make them 
different from equations whose solutions are globally meromorphic, such as the KdV equation. 
Second, we blend the techniques of algebraic geometry and weak solutions of PDE's to gain 
further insight into, and explicit formulas for, piecewise-smooth finite-gap solutions. 

The basic technique used to achieve these aims is rather different from earlier papers dealing 
with peaked solutions. First, profiles of the finite-gap piecewise smooth solutions are linked 
to certain finite dimensional billiard dynamical systems and ellipsoidal billiards. Second, after 
reducing the solution of certain finite dimensional Hamiltonian systems on Riemann surfaces 
to the solution of a nonstandard Jacobi inversion problem, this is resolved by introducing new 
parametrizations. 

Amongst other natural consequences of the algebraic-geometric approach, we find finite 
dimensional integrable Hamiltonian dynamical systems describing the motion of peaks in the 
finite-gap as well as the limiting (soliton) cases, and solve them exactly. The dynamics of the 
peaks is also obtained by using Jacobi inversion problems. Finally, we relate our method to the 
shock wave approach for weak solutions of wave equations by determining jump conditions at 
the peak location. 
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1 Introduction 



An important feature of many integrable nonlinear evolution equations is the nature of their soliton 
solutions. There are many examples of such solutions found in a variety of physical applications, 
such as nonlinear optics and water wave equations. Nonsmooth soliton solutions of integrable 
equations are now well known, and include solutions of the shallow water equation (|SWD with 
peaks, the points at which their spatial derivative changes sign (see Camassa and Holm [1993] and 
Camassa, Holm and Hyman [1994]). It was noted in Alber et al. [1994, 1995, 1999] that the spatial 
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structure of these "peakon" and finite-gap piecewise smooth weak solutions are closely related to 
finite dimensional integrable billiard systems. 

The purpose of this paper is twofold. First, we exhibit some of the special features of integrable 
PDE's that admit piecewise smooth weak solutions, which make them different from equations 
whose solutions are globally meromorphic, such as the KdV equation. Second, we blend the tech- 
niques of algebraic geometry and weak solutions of PDE's to gain further insight into finite-gap 
solutions that are piecewise smooth. 

Some history. Camassa and Holm [1993] described classes of n-peakon solutions for an integrable 
equation in the context of a model for shallow water theory. This work (see also Camassa, Holm 
and Hyman [1994]) contains many other facts about these equations as well, such as a Hamiltonian 
derivation of the equation, the associated linear isospectral eigenvalue problem and its discrete 
spectrum corresponding to the peakons, a steepening lemma important for understanding how 
solutions lose regularity, numerical stability, etc. Of particular interest to us is their description 
of the dynamics of the peakons in terms of a finite-dimensional completely integrable Hamiltonian 
system. In other words, each peakon solution can be associated with a mechanical system of 
moving particles. Calogero [1995] and Calogero and Francoise [1996] further extended the class of 
mechanical systems of this type. 

It is well-known (see, for example, Ablowitz and Segur [1981]), that solitons and quasi-periodic 
solutions of most classical integrable equations can be obtained by using the inverse scattering 
transform (1ST) method. This is done by establishing a connection with an isospectral eigenvalue 
problem for an associated operator that is often a Schrodinger operator. In some cases it involves 
a potential in the form of an entire function of the spectral parameter. Such an operator is called 
an energy- dependent Schrodinger operator. The scattering problem for the operators of this type 
was studied by Jaulent [1972] and Jaulent and Jean [1976]. 

On the other hand, in connection with certain A^-component systems of integrable evolution 
equations, Antonowicz and Fordy [1989] investigated certain energy dependent scalar Schrodinger 
operators. Using this formalism, they obtained multi-Hamiltonian structures for this class of sys- 
tems. 

Later, Alber et al. [1994, 1995, 1999] showed that in case of certain potentials, a limiting 
procedure can be applied to generic solutions, which results in solutions with peaks. The latter 
were related to finite dimensional integrable dynamical systems with reflections and were termed 
piecewise- smooth solutions, a terminology that hereafter we will adopt. This relation provides an 
efficient route to the study of finite-gap and piecewise soliton solutions of nonlinear PDE's. The 
approach is based on studying finite dimensional Hamiltonian systems on certain Riemann surfaces 
and can be used for a number of equations including the shallow water equation, the Dym type 
equation, as well as certain A^-component systems and equations in their hierarchies. 

Finite-gap solutions of the Dym equation were studied in Dmitrieva [1993a] and Novikov [1999] 
by making use of a connection with the KdV equation and with the aid of additional phase functions. 
Soliton solutions of Dym type equations were studied in Dmitrieva [1993b]. Periodic solutions of 
the shallow water equation were discussed in McKean and Constantin [1999]. The papers by Beals 
et al. [1998, 1999, 2000] used Stieltjes' theorem on continued fractions and the classical moment 
problem for studying multi-peakon solutions of the ( pW|) equation. Multi-peakon solutions have 
also been derived in Camassa [2000] by Gram-Schmidt orthogonalization. 

The main results of this paper. While our techniques are rather general and can be applied to 
large classes of A^-component systems, we shall illustrate them in detail for two specific integrable 
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PDE's. One of these equations is a member of the Dym hierarchy that has been studied by, amongst 
others, Kruskal [1975], Cewen [1990], Hunter and Zheng [1994] and Alber et al. [1995, 1999]. Using 
subscript notation for partial derivatives, this equation is 

U^xt + 2U^U^^ + UU^^^ - 2kU^ = . (HD) 

The other equation, derived from the Euler equations of hydrodynamics in a shallow water frame- 
work by Camassa and Holm [1993], is 

Ut + sua, = U^^t + U^x + UU^^^ - 2kU, . (SW) 

In both equations, the dependent variable U{x, t) may be interpreted as a horizontal fluid velocity 
and K is a parameter. 

Under appropriate boundary conditions, applying the limit k — > to ( |SW| ) leads to an equation 



that has peaked solutions. For equation ( HD ), such solutions exist also for k 7^ (for example 
periodic and finite-gap peaked solutions). 

By using the method of generating equations for nonlinear integrable PDE's, we reduce the 
equations to a Jacobi inversion problem associated with hyperelliptic curves. The solutions U{x, t) 
themselves are given by trace formulae, i.e., sums of coordinates of points on such curves. 

An important feature is that the corresponding Abel-Jacobi mapping is not a standard one. 
First of all, the holomorphic differentials that are involved do not form a complete set of such 
differentials on a hyperelliptic curve. Second, it involves a meromorphic differential. As a result, 
the image of the mapping turns out to be a non-Abelian subvariety (a stratum) of a generalized 



Jacobian. This also implies that the x- and t-flows of ( HD ) and ( ^W| ) are essentially nonlinear, 
i.e., they are not translationally invariant. Seen from the viewpoint of algebraic geometry, these 
nonstandard aspects constitute the main difference between shallow water and Dym type equations, 
and equations of KdV type and more generally equations from the whole KP hierarchy which lead 
to standard Abel-Jacobi mappings. 

The basic technique of the present paper is rather different from earlier papers dealing with 
peaked solutions. First, profiles of the finite-gap piecewise-smooth solutions are linked to certain 
finite dimensional billiard dynamical systems and ellipsoidal billiards. Second, after reducing the 
solution of the finite dimensional Hamiltonian systems on Riemann surfaces to the solution of a 
nonstandard Jacobi inversion problem, it is resolved by introducing new parametrizations. 

The philosophy that "justifies" procedures of this sort is that, in the end, by using the trace 



formulae, we obtain weak solutions of the PDE's ^HL\) and ( SW ) in the spacetime sense. This is 



regarded as equivalent to the validity of Hamilton's principle for these PDE's and is taken as a 
fundamental criterion for the definition of their solutions. It is worth emphasizing that Hamilton's 
principle naturally leads to weak solutions in the spacetime sense (and not in the spatial sense 
alone). We might also remark that even for billiards, one has to be careful about the sense in 
which solutions are interpreted. In the case of a point particle bouncing off a wall, for example, 
the equations of motion themselves do not rigorously make sense at the collision; what does make 
sense is the fundamental principle of Hamilton. This point of view of course is not new — see, e.g., 
Young [1969] and Kane et al [1999]. 



The contents of the paper. In §2, basic trace formulae and ^u-variable representations are used 
to establish a connection between solutions of the nonlinear equations and finite dimensional Hamil- 
tonian systems on Riemann surfaces. These representations describe finite-gap and soliton type 
solutions, as well as mixed soliton-finite-gap solutions. Then, solving the Hamiltonian systems is 
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reduced to Jacobi inversion problems with meromorphic differentials. These inversion problems are 
solved by introducing a new parameterization that yields a Hamiltonian flow on a nonlinear subva- 
riety of the Jacobi variety. The approach of recurrence chains used in this section is demonstrated 
in detail in the case of Dym-type equations. 

In §3 the geodesic motion and motion in the field of a Hooke potential on an ellipsoid Q are 
linked, at any fixed time t, to finite-gap solutions of ( |HD| ) and ( SW] ) equations respectively through 
trace formulae. In §4 it is shown how peaked finite-gap solutions of ( HP ) and ( SW ) equations 
arise in the particular limit of smooth solutions. Based on this, a connection to ellipsoidal and 
hyperbolic billiards is used to construct the peak solutions of equations ( HP ) and ( SW ) in the form 
of an infinite sequence of pieces, corresponding to the segments between impacts, glued together 
along peaks. The motion between impacts in the billiard problems is made linear on generalized 
Jacobians of hyper elliptic curves. 

By solving the corresponding generalized Jacobi inversion problem, we find theta-function so- 
lutions to the billiards, which thereby enables us to describe explicit peak solutions for the above 
PPE. We then extend the analysis from fixed-time peak solutions to time-dependent ones and show 
that the latter are described by an infinite number of meromorphic pieces in x and t that are glued 
along peak lines (surfaces) where the solution has discontinuous derivatives in the dependent vari- 
ables. We give theta-function expressions for the pieces and the peak surfaces. These formulae may 
be useful for stability analysis as well as for numerical investigations of the perturbed nonlinear 
PPE's. 

In §5 the Hamiltonian structure for the motion of the peaks of the finite-gap piecewise-solutions 
is obtained by using algebraic- geometric methods. Lastly, in §6 we relate our method to the shock 
wave approach for weak solutions of wave equations by determining jump conditions at the shock 
location. 



2 Finite-gap solutions 

In this section we will show that even on the level of finite-gap solutions, there are crucial differences 



between the KdV equation case and equations ( |IIP| ) or (SW). The same method can be applied to 
other equations forming the HP and SW hierarchy as well as to A^-component systems of nonlinear 
evolution equations which have associated with them energy dependent Schrodinger operators (see 
Alber et al. [1997]) . 

We will start by describing the algebraic geometrical structure of finite-gap solutions of equations 
(HP) and ( SWj ) related to a hyperelliptic curve of genus n, also called n-gap solutions. The same 



method can be applied also to the other equations forming the HP and SW hierarchy. 

For the HP equation such solutions were obtained in terms of theta-functions by Pmitrieva 



[1993a] (see also Pmitrieva [1993b]) and Novikov [1999]. For equation (SW) on a circle, the problem 
was discussed in Constantin and McKean [1999]. 

Lax pairs and recurrence chains. We now use the recurrence chain approach to develop a 



basic trace formula which establishes a connection between solutions of equation ( HP ) and finite 
dimensional Hamiltonian systems on Riemann surfaces, written in the so-called /i-variables rep- 
resentation. This representation describes finite-gap solutions, as well as their limiting forms of 
soliton-type. This representation also yields the existence of peakons in a special limiting case. For 
definiteness, we concentrate here on equation (|HP| ). Analogous results are available in the case of 
equation ( ^ ) (for details see Alber et al. [1994, 1995]). 
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The hierarchy of Dym equations is obtained from the Lax equations 
where the potential V [E ,x^tn) is written in terms of a complex parameter E in the form 

y{x.tn.E) = ^^^^^ (2.1) 

for a function M(x, t„) to be determined below. Assuming [L, to be a scalar operator, we choose 
An = Bndx — \B'n ^r some function Bn{E, x,tn) and obtain the following sequence of equations 
for V, 

Now we choose Bn to be a polynomial in E of degree n: 

n n 

Bn{x,t,E)=bo lliE - fik{x,t)) = J2bn-k{x,t) E'' . (2.3) 

fc=l k=0 



Substituting the expressions ( pHp and ( |2.3| ) into the generating equation (2^) and equating like 
powers of E, we obtain a recurrence chain for coefficients of B{x, t) which yields n-th equation of the 
Dym hierarchy. For example, putting ti = t and choosing n = 1, Bi{x,t,E) = bo{x,t)E + bi{x,t) 
yields the following chain 



1 

E^ : -6o =0 



E° : -b" + 26'oM + boM' = 
E-^ : 2b'^M + biM' = (2.4) 



After setting bo = 1 and using ( |2.1| ), we get 

M' = b'i' 

2b[M + biM' = (2.5) 

The first equation defines bi in terms of M, M = b'^ + k, with k a constant. Renaming hi = —U, 
so that 



M = + K, (2.6) 



and putting this into the second equation of the set (|2.5D results in equation (|HDD . (For further 
details about the hierarchies of ( pDD and ( pW] ), see for example Alber et al. [1994, 1995, 1999].) 
The method of generating equations is due to S. Alber [1979] and another exposition of it can be 
found in Alber et al. [1985, 1997]. 

We call ( j2.2| ) the "dynamical generating equations," because it generates a hierarchy of equations 
governing the dynamics of the dependent variable M{x,t). 



Remark. The flows where Bn is a polynomial E, as in the definition (|2.3| ) and in the example 
above, will in general lead to nonlocal equations, i.e., the evolution equation for M involves terms 
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that depend on nonlocal operators acting on combinations of M and its derivatives. This can be 



seen, for instance, in equation ( ^.51) where both bi and 6'^ require inverting (|2.6D to write U in 
terms of M. Thus, flows generated by polynomials Bn in E should properly classified as integro- 
differential evolution equations, rather than PDE's. In contrast, the choice of polynomials 'm.l/E 
for Bn leads to flows that are local, i.e., Mt only depends on combinations M and its (spatial) 
derivatives, and these flows are proper PDE's. This feature of equations of Dym (p3D|) or shallow 



water ( ^W| ) type is somewhat different from other completely integrable PDE's like the KdV or 
Sine-Gordon equation. Equations ( |HD| ) and ( |SW| ) possess "open ended" hierarchies: the recurrence 
chain can be extended from negative to positive powers of E, by choosing Bn in ( ^.2[) to be a rational 
function of the parameter E. The case when the chain includes only negative powers of E is in fact 
the one most studied in the literature (see, e.g., Dimitrieva [1993a], Novikov [1999] for the case of 
Dym equation). 

Now let us consider the stationary flow for the n-th equation of the hierarchy, which is obtained 
by dropping the time derivative of V in the left-hand side of (|2.2|) . By definition stationary equation 
describes a finite-dimensional system for the coefficients of Bn and is equivalent to the 2x2 Lax 
pair 



^Wn{E) = -[Wn{E),C{E)l or 



^ + C{E),Wn{E) 



0, (2.7) 



Wn{E) — ( I „2 ^ .J 1 R/ ) ' ^ — { M n 
-^Bn + Bn^ ^BnJ \^ U 

The matrix WniE) undergoes an isospectral deformation. Hence the spectral curve 

T = {\Wn{E) - zl\ = 0} 

is an invariant of the stationary flow. The curve is hyperelliptic and can be represented in the form 

r = {w' = fiCifi)} (2.8) 

where z = wE and 

C{E) = E {-B'nBn + \B'f^ + BlM. (2.9) 
Since Bn is a polynomial of degree n, C{E) becomes a polynomial of degree (at most) 2n: 

2n 2n 

C{E) = Cj E^ = C2n J{{E- mfc), (2.10) 
j=0 fc=0 

for some constants m^. A; = 1, . . . , 2n. In this case the curve V has genus n and we set coefficient 
C2n to be a negative number: C2n = —Lq. We shall refer to ( ^ ) as the stationary generating 
equation. 

Equating like-powers of E in both sides of the stationary generating equation yields first integrals 

^'":C2„ =-b'i+M, 
^2-1 : C2n-i= - hb'i - 62 + ^{b[f + 26iM, ^2 



E^ : 

^^0 



E° : Cn =2blM. 
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Let us consider divisor of points Pi = wi), . . . , P„ = {fin,Wn) on T. Substituting ( ^^ ) into 
(|2.9|) and setting E = ixi, . . . ^ successively, one gets the following system of equations describing 
evolution of the points under the stationary flow, 



dm _ 



(2.12) 



where 



2n 



Ulr 



(2.13) 



r=l 



In the case of equation ( pW| ), this should be replaced by 

2n+l 



rrir 



(2.14) 



r=l 



We now proceed to describe finite-gap solutions of equation ( |HD| ) and the other equations from 
its hierarchy. According to a general theory (see, e.g., Dubrovin [1981], Belokolos et al. [1994], for 
any fixed t, the x-profile of an n-gap solution of an integrable PDE satisfies the n-th stationary 
equation of the hierarchy. Hence, n-gap solutions M{x,tk) of k-th. equation of HD hierarchy must 



satisfy the stationary generating equation ( |2.9| ) represented by the Lax pair (2.7), as well as the 
dynamical generating equation 



dV 

dtk 



I d^Bk ^ 2^^^V I B — 
2 dx^ dx ^ dx 



V 



M{x,tk) 
2E ' 



(2.15) 



where the coefhcients of Bk{E) are found recursively. Notice that the latter equation is equivalent 
to the matrix commutativity relation 



ox atk 



0, 



where 



Wk{E) 



Bk 

M 1^, ), 



and C is defined in ( p. 7]) . The compatibility of conditions ( 2.16| ), and 
Lax pair 



_d_ 

dtk 



Wn{E) = -[Wn{E),Wkm. ^eN, fc^n. 



(2.16) 
(2.17) 

leads to the following 
(2.18) 



For k = n, we replace ( p. 18 ) with the Lax pair (2.7) thus identifying t„, with x. 

The (l,2)-entry of the matrix equation ( ^.18 ) implies the following t^-evolution of the polynomial 

Bn{E) 



dBr, 

dtk 



dBn dBk 



dx 



dx 



k ^ n. 



(2.19) 
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In case k = n this relation is replaced by 



dB 



vbo 



dB 
dx 



where f is a constant, which can always be eliminated by rescaling 

Expanding the right hand side of ( 2.ig| ) in E and using the condition that it must be a polyno- 
mial of degree n — 1, we find 



BkiE) 



1 



En^k 



-Bn{E) 



(2.20) 



where [ ]+ denotes the polynomial part of expression. 

As follows from the first equation in (|2.11| ), M = C2n + On the other hand, according to 
formula ([2.3|), hi = — X^^Li Finally, using the definition ( [2.6[ ) of M in terms of the solution U 
and integrating twice with respect to x, we obtain 



U 



n ^ 

X] + 9 ('^ - '^2n,) x"^ + KiX + K2 



(2.21) 



where Ki and K2 are constants of integration. If we assume that all the variables /Xj are bounded, 
which is related to the choice of sign of the leading order coefficient C2n, then bi is a bounded 
function of x. To find bounded solutions U{x, t) of the PDE, we set 



G 



2n 



and Ki = 0. 



Hence, when the above requirements are imposed, we see that the leading order coefficient of the 
polynomial C{E) must coincide with the parameter k of the PDE. 

The Dym equation ( |IID| ) is invariant under the Galilean transformation 

x = x + K2t, i = t, U = U - K2, 



so that the constant K2 can always be eliminated from expression ( 2.21| ). Therefore, under the 
boundedness conditions above, and up to a Galilean transformation, we assume that the finite-gap 



and soliton solutions of the Dym equation ( IHDD is reconstructed in terms of the root variables n's 
by the "trace" formula which in case of equations ( |IID| ) and ( SW ) have the form 



U{x,t) = ^fJ^i 



m . 



(2.22) 



i=l 



Here m is a constant, which equals zero in the case of equation 

Through ( 2.22| ) a solution of the system ( |2.12| ) allows to construct the instantaneous profile of 
U{x, •) from a set of initial conditions fj.i{x, •) = ^i(0, •) G [m2i,m2i+i], i = 1, . . . ,n. Here the "dot" 
notation stresses the fact that time t is just a parameter in this system. 

, successively, and 



On the other hand, substitution of (|2.3D into ( |2.19[) , setting E = fii, 
taking into account expressions ( 2.12| ) results in the following ife-evolution equations for jii 



dtk 



(2.23) 
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where, in view of (2.3) and (2.20), for k = 1, 

Bkim) = ReSs=o ■ 



,n 



1 



i.e., up to the sign, the k-th. elementary symmetric function of {/ii 



■ A*n} \ In the case A: = 1, 



1, 



, n, 



S = /Xl h /4,i 



(2.24) 



solution of which produces /t's, and hence the PDE's solution U , at any (later) time t. We notice 
that for k > n, the derivatives d/dtt are linear combinations of d/dti, . . . , d/dtn- 

Expressions ( 2.12 ), ( 2.23| ), and ( [2.12 ) provide the so-called ^-variables representation for the 
finite-gap solutions of an evolution equation. They are the analogs of Drach-Dubrovin equations 
which describe evolution of points on the spectral hyperelliptic curve in the case of the KdV 
equation. (For further details see Dubrovin [1975], Drach [1919], Alber et al. [1994, 1995, 1999], 
Gesztesy et al. [1996], and Alber and Fedorov [2001].) 

With the initial conditions chosen the right-hand-side of system ( [2. 121) is real, and the derivative 
of /tj changes sign when /ij reaches the end points of its gap, /ij = m2i or fii = m2i+i, corresponding 
to change of the sheet of the spectral curve F. Thus each variable /z undergoes (real) oscillations 
between the end points of a gap (so that the resulting PDF solution U{x,t) remains real). 

Remark. The condition that the root variables /x's are real (or, equivalently, their initial conditions 
are chosen as described above), while certainly sufficient to assure reality of the PDF's solution U 
resulting from (2.21), is clearly not necessary (namely, some of the /i's could occur in conjugate 
pairs). A wider class of real solutions U could be constructed by relaxing the reality assumption on 
the /i-variables. However, a thorough discussion of the reality condition for U and its implications 
for the root variables, while certainly desirable, lies beyond the scopes of the present paper, and it 
will be addressed in future work. 

By rearranging and summing up equations ( |2.12| ) and ( |2.24| ), ( p. 23 ), one obtains the following 
nonstandard Abel-Jacobi equations 

fJ-idfii _ [ dtk 
dx 



n 



k 
k 



n, 



(2.25) 



which contain (n — 1) holomorphic differentials and one meromorphic differential on F. Thus, the 
number of holomorphic differentials is less than genus of the Riemann surface, which implies that 
the corresponding inversion problem cannot be solved in terms of meromorphic functions of x and 
ti, . . . ,tn-i (see e.g., Markushevich [1977]). 



Finite-gap stationary flows in x. Let us first consider the x-fiow by fixing time variables in 
( 2.25| ): tk = =const, fc = 1, . . . , n — 1, so that dt^ = 0. Now introduce a new spatial variable xi 
defined as follows 

— Hi - ■ ■ Undxi. (2.26) 
^0 

In view of the well-known Jacobi identities 

' l/(m 

//• _ J ^ ^2.27) 



1 



fin) k = -l, 

k = 0, . . . ,n 

k = n — 1, 
k = n, 
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the equations ( |2.12| ) give rise to the fohowing system 

r + 0^ A; = 1, 



E 



MO 



yjR{li) \ (t^k k = 2,...,n, 



(2.28) 



where (pi, . . . ,(j)n are constant phases which depend on as on parameters. 



Equations (2.28) include n holomorphic differentials on T and determine the standard Abel- 



Jacobi map of the symmetric product F^") of n copies of T to the Jacobi variety (Jacobian) Jac(r). 
Thus, the flow generated by the system ( |2.12| ) is made linear on Jac(r) after introducing the 
reparametrization ( ^.261) . By using standard methods (see e.g., Dubrovin [1981] or Mumford [1983]), 
the map can be inverted, resulting in expressions for algebraic symmetric functions of /i-variables 
in terms of theta-functions of n arguments which depend linearly on xi and, in a transcendental 
way, on as parameters. Then, by using the trace formula ( ^321) , one obtains a theta-functional 
expression for [/ as a function of xi,t^, U = U{xi,t^). 

On the other hand, substituting the theta-functional expression for the product • • • /in into 
( 2.26| ) yields a quadrature. By solving it, one finds x as a meromorphic function of xi which depends 



on to as a parameter. However, the inverse function xi{x,to) is no longer meromorphic in x. 

Finally, the composition function C/(x,to) = U{xi{x,to),t^) gives a profile of the finite-gap 
solutions of the ( |HD| ) or ( ^Wj) equation (for explicit theta-functional expressions U{xi,tQ), x{xi,to) 
see Alber and Fedorov [2001]). Notice that as seen from ( |2.26| ) and ( |2.28| ), the original x-flow is 
also made linear on Jac(F). However the straight line motion is not uniform. 

The transformation ( |2.26| ) involving x and xi coincides with a change of variable in the well- 
known Liouville transformation (see, e.g., Verhulst [1996]). 

Finite-gap flows in tk- Now let us fix the coordinate x = xq as well as all the times ti, . . . , tn-i 
but tk- Then introduce a new time variable i defined by 

dU = T^T^.di, (2.29) 
where are the elementary symmetric functions oi ni, . . . , fin such that 

(S - /Zi) • • • (S - fin) = + + • • • + S^S,. 

Applying again the identities ( ^^ , from (|]2|) and (|2^ ) we arrive at the following canonical 



Abel-Jacobi mapping 

r"^ fi'-' dfi _ f = i + s = i, 



2y/R{fl) { ips = Ss,ktk + 4>s s = 2, ...,n. 



(2.30) 



where (pi, . . . ,(j)n are constant phases which depend on xq and the rest of times t; as on parameters, 
and 5ij is the Kronecker delta. 



As a result of inversion of ( 2.30| ), elementary symmetric functions of /i's and therefore the 



solution of equations (HD) and ( SW) ) can be found in terms of theta-functions of n arguments 



which depend linearly on ipg. This means that the arguments depend linearly on t, as well as on 
the original time tk- However, t itself depends on tk in a nonlinear way. 

Indeed, to describe the relation between t and t^, we substitute the theta-functional expressions 
for the symmetric functions S„ = //i • • • /tn and into ( 2.29| ). As a result, in contrast to the 



quadrature (2.26) relating x and xi, we now get a differential equation of the form 



^ = F{tk,i\xQ), 
dt 
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where is a transcendental function of t, t and the parameter xq. It can be shown that the equation 
involves a transcendental integral. 

Remarks. 

1. In contrast to the xi- and x-flows considered above, the flows generated by ( 2.23| ) (t^-flows) 



including (|]2|), are nonlinear flows on the Jacobi variety Jac(r). From the point of view of 
algebraic geometry, this phenomenon constitutes the main difference between solutions of such well 



known equations as KdV and sine Gordon equations and equations of (HD) or (SW) type. 

2. The problem of inversion of the full nonstandard Abel mapping defined by ( ^ ) can be also 
studied by using a generalized Jacobian of the curve F. Namely, one has to extend the mapping by 
including an extra holomorphic differential on F to get a complete set of such differentials. As a 
result of this procedure, one gets a flow on nonlinear subvarieties (strata) of generalized Jacobians. 
The complete algebraic geometrical description and explicit formulae are presented in Alber and 
Fedorov [2001]. 

3 Flows on n-dimensional quadrics and stationary n-gap solutions 
of the ( [HDl) and ( |SW|) equations. 

Consider a family of confocal quadrics in W^^^ = {Xi, . . . , Xn+i) 

( X"^ X"^ 1 

Qis) = \ — + ■■■ + ^±^ = lL sGM, < a„+i < ai < ••• < a„. (3.1) 

[ ai - s a„+i - s J 

The elliptic coordinates /^i, . . . ,/in+i can be defined in M""'"^ in a standard way (see, e.g., Jacobi 
[1884a]) as follows. The condition s = c determines the quadric Q{c) on which one of the coor- 
dinates, say Hn+i, equals c, and the other coordinates fii, . . . , fin are elliptic coordinates on Q{c) 
defined by relations 

X| = (a,-c)-iP^p^, , = l,...,n + l. (3.2) 

In the sequel without loss of generality we assume c = 0. 

It is well-known that the problem of geodesies on the ellipsoid Q = Q{0) is completely inte- 
grable (Jacobi [1884 a,b]). Moreover, as noticed by Jacobi himself and later by many other authors 
(see e.g. Ranch- Wo jciechowski [1995]), there exists an infinite sequence of integrable generaliza- 
tions of the problem describing a motion on Q in the force field of certain polynomial potentials 
Vp{Xi, . . . ,Xn+i), p G N of degree 2p. The simplest integrable potential is the quadratic Hooke 
potential or the potential of an elastic string joining the center of the ellipsoid Q to the point mass 
on it: 

Vi = |(X2 + ... + z2+i), a = const. 

In this case in terms of the ellipsoidal coordinates, the total energy (Hamiltonian) takes the Stackel 
form: 

8^ $(^,) [dx) +2 1./^^ + --^*' 



12 



where 

= - ai) • • • (/i - a„+i) 

and X denotes time. After fixing constants of motion, the system is reduced to the Abel-Jacobi 
equations 



^ /^FTT — V ^inX + (f>i, i = l,...,n, (3.3) 
TZ{fi) = -fi^{fi)[co{n - ci) ■ ■ ■ (/i - Cn-i) - fJ^"], Co, . . . ,c„_i = const, 

where (pi, . . . , (j)n are constant phases and ci, . . . , Cn-i are constants of motion. 

Notice that for a = the order of the polynomial TZ{fi) is 2n + 1, whereas for o" 7^ it is 
2n + 2. The case o" = corresponds to the free (geodesic) motion on Q. cq is the constant in 
the first integral {X,X) and the remaining constants admit a clear geometric interpretation: the 
tangent line to a geodesic is also tangent to the fixed confocal quadrics Q{ci), . . . , Q{cn-i) (Chasles 
theorem) . 

Now notice that equations ( |3.3| ) are equivalent to the system ( 2.25| ) with dt = describing 



stationary ( HP ) and (|SW| ) equations, provided we identify the roots of the polynomial TZ{fj,) with 



those of the odd order polynomial ( 2.13| ) (for o" = and Lq = 1) and of the even order polynomial 



( 2.14| ) (for fj = 1) respectively. The equivalence also holds when some of the parameters Oj in ( ^.3| 



are negative, which correspond to the motion on a hyperboloid. For concreteness, we shall consider 
only the case of ellipsoids. Taking into account the trace formula ( ^.22| ), we arrive at the following 
theorem 

Theorem 3.1 The geodesic motion and motion in the field of a Hooke potential on the ellipsoid 
Q are linked, at any fixed time t, to the n-gap solutions of ( \HL\ ) and ( ISU]) equations respectively 
through the trace formula \2.23^ ). Namely, if the roots of the polynomials R{ij,) in ( 2.13^ ) or ( [g. I4J 



coincide with the roots of TZ{fJ,) in ( \3. jj ), the profiles of such solutions are given by the sum of the 



elliptic coordinates of the moving point on Q with addition of (—xn) in case of equation (S]^') 



For the geodesic flow on Q {a = 0) and equation ( |HD| ), this result was obtained in Alber and Alber 
[1985], Cewen [1990], and Alber et al. [1995]). 

As with equation ( 2.25 ), under the change of parameter ( p. 26 ), equations (^^) reduce to those 



containing holomorphic differentials only and having the same structure as ( |2.28|) . By inverting 
the corresponding Abel-Jacobi mapping, one obtains explicit expressions for elementary symmet- 
ric functions of /ij and, in view of ( |3.2| ), for the Cartesian coordinates Xi, . . . ,Xn+i in terms of 
theta-functions of the new parameter xi (for the case of the geodesic flow, see Weierstrass [1844], 
Moser[1978], and Knorrer [1982]). In the case n = 2, the change of parameter ( [2.26 ) was first 



applied by Weierstrass [1844] to solve the classical Jacobi geodesic problem on a triaxial ellipsoid 
(Jacobi [1884a, 1884b]). 



Billiard dynamical systems and piecewise-smooth weak solu- 
tions of PDE's 



In this section it is first shown how peaked finite-gap solutions of ( HDD and ( |SW| ) equations arise in 



the limit rrii 0, where mi is the smallest root of the polynomial R{E) in equations ( |2.12| )- (|2.24| ). 
Then a connection to ellipsoidal and hyperbolic billiards is established. 
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Ellipsoidal billiards and generalized Jacobians. Suppose that one of the semi-axes of the 
elUpsoid Q tends to zero, namely, a^+i — > 0. In the limit, Q passes into the interior of (n — 1)- 
dimensional ellipsoid 

Q = {X^/ai + ■■■ + Xl/an = 1} e W", = {X,,. . . ,X„). 

The elliptic coordinates ni, fin on Q transform to elliptic coordinates in M" giving 

nkl^^p^, j = l,...,n, (4.1) 

which appear as the corresponding limits of (p.2|). 

Then the motion on Q gets transformed into billiard motion inside the ellipsoid Q. Geodesies 
on Q pass into straight line segments inside Q, whereas the points of intersection of the geodesies 
with the plane {Xn+i = 0} are mapped into impact points on Q with elastic reflection. Also, the 
motion on Q under the Hooke force passes to the motion inside Q under the action of the Hooke 
force with the potential V = <7{Xf + • • • + X^)/2. However, in contrast to cases a = or o" < 0, 
for fj > (an attracting Hooke potential), for the trajectory to reach Q the total energy h must be 
sufficiently large. Namely, there ought to exist a positive e such that inside Q the following double 
inequality holds 

h + a{Xf + ■■■ + Xl)/2 > e > 0. 

Under this condition, the motion on Q transforms to billiard motion inside the ellipsoid Q again 
having impacts and elastic reflections along Q. Thus, we have "an ellipsoidal billiard with the 
Hooke potential" which is described by the mapping B : (x, v) (x, v), where x, v E M" are the 
Cartesian coordinates of a point on Q and the starting velocity vector respectively, while (x, v) are 
the coordinates and the starting velocity at the next impact point. Following Fedorov [2001], the 
mapping has the form 

X = — i[(cr-(v,a-V))x + 2(x,a-V)v], 

V = — [(cj - (v,a-V))v - 2cj(x,a"V)x] + ga^^i (4.2) 

V 

= [(a — (v, a~"'^v))(v + ga~^:x.) + 2(x, a~^v)(ga~^^r — ax)], 

V 

^= VMx,a-iv)2 + (a-(v,a-iv))2, ^ = ^^1^;^. 

(x, a "^x) 

Notice that in the limit <t — > this reduces to a standard billiard mapping given in Veselov [1988] 



2(x,a-^v) _ 2(v,a~^i) 

x = x T^^-, v = v+— KTT^a X. 

(v,o ^v) (x, a ^xj 

The mapping ( ^.21 ), as well as the billiard limits of the motion on Q with the higher order potentials 
Vp(Xi, . . . , X„+i) (X„+i = 0) are completely integrable. 

In the limit a^+i — > and after using the change of variable ( p. 26 ), the Abel-Jacobi equations 
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|3.3| ) are transformed as follows 



k=l 



(pi = const, i = 1, . . . ,n — 1, 

} (4.3) 



Pil^) = -(^ - ai) • • • (/i - an) [coifj. - ci) - c„_i) - cr^"] 

This system contains n — 1 holomorphic differentials on the Riemann surface C = {w'^ = p{n)} of 
genus g = n — 1 and one differential of the third kind having a pair of simple poles Q_ , Q+ on 
C with /i(Q±) = 0. Here again (pi, . . . , (pn are constant phases and cq, . . . , c„_i are constants of 
motion. The elliptic coordinates /Ui, . . . , represent the divisor of n points Pi = {fii,Wi) on C. 

The equations ( |4.3| ) describe a well defined mapping of the symmetric product C^^"*"^^ to 
Jac(C, Q-, Q+), the {g + l)-dimensional generalized Jacobian of the curve C with two distinguished 
points Q± . The later is obtained from the genus n curve w = n{fi) in (O) as a result of confluence 



of two Weierstrass points (a„+i 0) and regularization: cutting out the double point and gluing 

Q-,Q+- 

The generalized Jacobian is a noncompact algebraic variety which is topologically equivalent to 
the product of the customary (7-dimensional Jacobian variety Jac(C) with complex angle coordinates 
(pi, . . . ,(pg and the cylinder C* = C\ {0} (for the definition and description of generalized Jacobians 
see, among others, Serre [1959], Previato [1985], Gavrilov [1999], and Fedorov [1999]). 

As follows from ( [4.3[ ), the geodesic and the potential billiard motion parameterized by xi is 
represented by a straight line flow on Jac(C, Q_, Q+), which is directed along the real section of 
C* and leaves the coordinates (p on Jac(C) invariant. 

As we shall see below, for the odd order and even order polynomial -R(/u) the solutions to the 
generalized inversion problem ([4. 3D have different structures. 



Solutions in terms of generalized thet a- functions. First we concentrate on straight line 
billiards corresponding to the case cr = when the curve C has one infinite point 00. Fix a canonical 
basis of cycles Ai, . . . , Ag, Bi,...,Bg on C and let (Di,...,u)g be the dual basis of normalized 
holomorphic differentials on C and zi, . . . ,Zg be corresponding coordinates on the universal covering 
of Jac(C). There exists a unique g x g constant normalizing matrix D such that 

Dkj A*"'"^ dfi 

^k = Z^ 7= — , Zk=yDkj(pj, k = l,...,g = n-l. (4.4) 

j^i Vp(.i^) J^i 

Let us also introduce a normalized differential of the third kind ^Iq having simple poles at Q± with 
residues ±1 respectively: 



^0 = r-TT + 5^ l^'^^k, \/p(0) = yjai • • • a„ • ci • • • c„,_i, (4.5) 
/^VP(^) k=i 

where f3k are unique constants such that r^o has zero A-periods on C. Then the last equation in 
|) can be represented in the following form 

^ / 9-Q = Z, Z = 2^/^Xl + const. (4.6) 

k=i ''t'o 
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Notice that in case of the ehipsoidal bihiards sj R{0) is always real and hence Z is also real. Let us 
also choose the base point {fiQ, wq) of the mapping (^^) to be an infinite point oo G C. According to 
Fedorov [1999], the solution of the problem of inversion ([1.3D together with ( |4.1| ) yields the following 
expressions for the Cartesian coordinates Xi of the point moving inside the ellipsoid Q: 



Z = {Zi,. . . ,Zn~lf , Z = 2yjR{Q)xi + Zq, Z, Zq = COUSt, 

+ _ rQ+ _ \ ^ 

wi, . . . , / Wo I G C^, Hi = const. 



CXD J OO 



These expressions involve quotients of generalized theta-functions, where 0[A + ?/(j-)](z) and 0[A](2;) 
are customary theta-functions associated with the Riemann surface C with appropriately chosen 
half-integer theta-characteristics r/j-j) (A is the half-integer theta-characteristic corresponding to 
the vector of Riemann's constants). The vector q coincides with the vector of i?-periods of the 
meromorphic differential Oq. The constant factors Kj depend on the parameters of the curve C 
only. (For the definition and properties of the generalized theta-functions see e.g., Belokolos et al. 
[1994], Gagnon et al [1992], Ercolani [1987], and Fedorov [1999].) 

The expressions (l4.7|) describe a straight line segment in M" (C") with z playing a role of a 
constant phase vector which defines the position of the segment. When one of the /i-variables, say 
/ii, equals zero, the corresponding point Pi = (^ui, ^/R{JM)) on the curve C coincides with one of 



the poles Q_, Q+ of the differential Oq. Then, as follows from the mapping (O) and (p~q), xi and 
Z become infinite. On the other hand, in view of ( [l.lD , at this moment the moving point in R" 
meets an ellipsoid Q. 

It follows that as xi and Z change from — oo to oo along the real axis, the expressions ( |4.7| ) have 
finite limits, giving the coordinates of two subsequent impact points on Q. Notice that Xi{oo,z) 
have the same values as cx), z -\- q). Hence the next segment of the billiard trajectory is given 



by (4/7) with z being replaced by z + g. This yields the following algebraic-geometrical description 



of the billiard motion (see also Fedorov [1999]). 

Theorem 4.1 As the point mass inside Q approaches the ellipsoid, the point Pi on C tends to the 
pole Q+. At the moment of impact, Pi jumps from back to Q_, whereas the phase vector z is 
increased by q. The process repeats itself for each impact. 

Using this property and by applying induction, from ( [4.7D the coordinates of the whole sequence 
of impact points are found in the form 

g[A + r/(,)](zo+iVg) 

emzo + Nq) ' ^ =!'•••'" (4.8) 

where G N is the number of impacts and the phase vector zq = {zio, . . . , Zgo)^ is the same for all 
the segments of the billiard trajectory. 

These expressions depend on customary theta-functions only and, as functions of zq, are mero- 
morphic on a covering of the Jacobian variety of C. They have also been obtained by Veselov [1988] 
by using a factorization of matrix polynomials (see also Moser and Veselov [1991]). The work of 
Veselov is closely related to the discretization of mechanics that preserves the integrable structure. 
The numerical implementation of Veselov's procedures was given in Wendlandt and Marsden [1997], 
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a discrete reduction procedure in Marsden, Pekarsky and Shkoller [1999], Bobenko and Suris [1999] 
and an extension to PDE's in Marsden, Patrick and Shkoller [1999]. 

The generalized Abel map ( |4.3[ ) yields expressions in terms of generalized theta-functions for 
the elementary symmetric functions of the variables In particular, following Fedorov [1999], one 
obtains 



/"I • • • /"n 



d,,dvioge[A]{z,z) 

e-^/29y0[A](z - q/2) + e^/'^dv0[A]{z + q/2) 



e-^/^e[A]{z - q/2) + e^/^e[A]{z + q/2) 



(4.9) 



where 



9[A]{z, Z) = e~^/^e[A]{z - q/2) + e^/^9[A]{z + q/2), 

dv 



Z = 2y^xi + Zo, 



Vi— + + 

OZi dZn 



(4.10) 



V = {Dig, 



and where V is the last column of the normalizing matrix D defined in 
The phases z and Zq are the same as in (|4.7| ). As follows from ( [4. 91) , for xi, Z — > ±00, the product 
fj,i - • • fin tends to zero, as expected. Taking the integral (2.26) with Lq = 1 yields 



x{xi,z) 



l-ii - ■ ■ fin dxi = dv log 6{z, Z) + const 



e-^/^dve[A]{z - q/2) + e^/^9yg[A](z + q/2) 
e-^lH{A\[z - q/2) + e^lH{A\{z + q/2) 



+ const. 



(4.11) 



It follows from this expression that the original parameter x has finite values as x\ — > ±00 and 
x{qo,z) has the same value as x(— oo,z + q). Now, substituting in ( |4.11| ) Z = —00, Z = 00, by 
induction, we find the length of the A^-th segment of the billiard trajectory in the form 



x{N) - x{N - I) 



dve[A]{zo + Nq) dve[A]{zo + Nq 



e{A]{zo+Nq) e[A]{zo+Nq-q) 



N £ 



(4.12) 



zq being the same as in 

As a result, the solution Xi(x), x G M, of the continuous geodesic billiard problem should be 
viewed as consisting of infinite number of pieces each parameterized by xi E (—00, 00) and given 
by (4.7) and ( 4.11 ). These pieces are obtained by iteratively adding vector q to the phase z in { \i.T\ j 
and (|4.lT) and they are glued together at the impact points corresponding to xi = ±00. 



Now we turn to the ellipsoidal billiard with the Hooke potential (a = 1). In this case the curve C 
appearing in ( [4. 3D has 2 infinite points at ±00. We again introduce normalized differentials u)fc, ^Iq, 
and coordinates z^, Z according to ( ^.4| ) and ( [4.61) . Let the base point of the mapping ( |4.3| ) be one 
of the Weierstrass points of C, say fiQ = a„. Then, instead of ( [4.7] ), the inversion of the generalized 
mapping (|4.3|) yields the following expressions for the squares of the Cartesian coordinates of the 
mass point moving inside an ellipsoid Q: 



Xfixuz) 



6i2[A + r?(,)](z,Z) 



1, . . . ,n. 



e[A]{z - q/2, Z - S/2) e[A]{z + q/2, Z + 5/2) 

z = {zi, . . . , Zn-i)^ , Z = 2y/ p{0)xi + Zq, z, Zq = const. 
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(CJI, . . .,OJg)^ 



[uJi, . . . ,u;g) , 



s 



OO-)- 



(4.13) 
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where 0[A](2;, Z) is defined in (|4.10D and 

e[A + r?(,)](z, Z) = e~^lH{l\ + 7?(,)](z - g/2) + e^/2^[A + + (//2) (4.14) 

Here k- are constants, and ^ p{0) is the same as in ( |4.5| ). Similarly to ( [4.7[ ), as xi and Z pass from 
— oo to oo, Xf{xi, z) tend to finite values resulting in the squares of the coordinates of subsequent 
impact points on Q. Thus, expressions (|1D describe a segment of trajectory of the billiard in the 



field of the Hooke potential between two impacts. After each impact the phase vector z changes 
according to Theorem 4.1. Then, by using induction, the sequence of impact points is described as 
follows 

''^^'^> '''e[/\]{zo-q + Nq)e[A]{z, + q + Nqy 

A^GN, zo = (2io,---,^go)^ = const. (4.15) 

Apparently, this theta-functional solution for the billiard with the Hooke potential was not previ- 
ously known. Lastly, we find the following expression for x 

x(xi, z) = const + log ^^^^^^ ~ ^ ~ ^J'^^ , Z = 27^X1 + const, (4.16) 
e[^\{z + q/2,Z + 3/2)' V/^W i 

which, for xi ±00 and Z ±00, has finite limits determining x for two subsequent impacts. 
Then, using the expression ( 4.1C| ), by induction, we express a x-interval between the impacts in 



terms of the customary theta-function: 

x(N) - x(N - 1) = log ^[A](^0-^/2 + iVg) _ Omz,-q/2 + Nq-q) _ 

x^iy) x^iV ij ^°g^[^](,^ + ^Y2 + Arg) '''^e[A]{zo + q/2 + Nq-q) ^^"''^ 

We emphasize that, in contrast to the geodesic billiard, for the billiard in the potential field the 
"time" X is not proportional to the length of a trajectory. 

Stationary finite-gap peaked solutions. Now we return to the finite-gap solutions of equations 
(HD) and (SW). Notice that under the limit mi — > the mapping (|]2|) takes the form (gj) with 
p{p) being a polynomial of degree 2n — 1 and 2n respectively. 
The trace formula ( p.22| ) and relations ( [4.1D yield 

n n 

C/ = ^ A| ^ + m. 

j=i i=i 

Then solution to the billiard problems ([4.7D-( ^rT7[ ) provide solutions U (x, to) for the above equations 
which consist of infinite sequences of smooth pieces each one corresponding to a segment between 
two impacts. The impacts themselves give peaks of C/(x,to)- This leads to the following theorem. 



Theorem 4.2 1). At any fixed time t = to, finite-gap peaked solution of the equation (HD) consists 



of an infinite number of pieces Ui~f{x,to), A^ G Z glued at peak points. Let p(/i) he any polynomial 
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with distinct roots ai, . . . , a„. Then, for any N, every piece is given by the following pair of theta- 
functional expressions parameterized by xi G M, 



n 



j=l i=l 

_ e-^l^dyOmz^ - g/2) + e^l^dye[^]{z^ + g/2) 
^ ^ " e-^/^emzM - q/2) + e^/20[A](z^ + q/2) + ^^"'^^ 
ZN = ZQ + Nq(^C^-^, Z = 2y^xi + Zo, 



where Xj{xi, z) and q are given by (4-V and zq, Zq, xq are constant phases of the solution depending 
on to, which are the same for any piece. The length of the N-th piece equals 

dve[^\{zo + Nq) _ dve[A]{zQ + Nq-q) 
e[^]{zo + Nq) em{zo + Nq-q) ' ^ ' > 



2). At any fixed time t = to finite-gap peaked solution to equation ( SU\ ) consists of an infinite 



number of pieces UN{x,to), N £ Tj which are glued at peak points. The pieces are given in the 
following parametric form 

n n 

Un = J2xf{xi,ZN) + J2a^ + m, zn = zo + NqeC'-\ (4.21) 
j=i i=i 

m{zN-q/2,Z-S/2) 



x(xi,to)=iog ,r:p y,: +^o, z=2yWu+zo, (4.22) 

9[A]{zN + q/2,Z + S/2) 

where Xj{xi,z) are given by (j^4-lV and zo,Zo,xo are constant phases which depend on to. The 
x-length of N-th piece equals 

emzo-q/2 + Nq) _ ^[A](.o - ^/2 + iV, _ ,) _ 
°^^[A](zo + g72 + iVg) ^e[A]{zo + q/2 + Nq-q) ^ ""^ 



When in the polynomials ( 2.13| ) or ( |2.14| ) mi = and 772-2 tends to zero, the distance between 



subsequent peaks of a profile tends to zero and in the limit the peaks coalesce. (Notice that this is 
done for a fixed t.) The solution U{x,to) for this limiting case is smooth. 

Remark. It is known (see, for instance, Fedorov [1999]) that there are special degenerate umbilic 
billiard solutions of the classical billiard problem (without a potential) that have straight line 
segments meeting n — 1 fixed focal conies of Q between any subsequent impacts and, as x ^ ±00, 
the billiard motion converges to simple oscillations along the largest axis of the ellipsoid. This 
corresponds to the confluence of the roots of the polynomial p{fi) in ( |4.3| ), 

Cl = ai, Cn-l = On-l- 

As a result, the hyperelliptic curve C becomes singular of arithmetic genus zero and the asymp- 
totic billiard motion is described in terms of tau-functions. The corresponding asymptotic peaked 
solutions of equations (HD) and (SW) are given in Alber and Fedorov [2001]. 
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Time-dependent piecewise-meromorphic solutions. Now we pass to global algebraic geo- 
metrical description of the finite-gap peaked solutions. After setting mi — > 0, the system (|2.25| ) is 
formally reduced to the following Abel-Jacobi mapping 



1 dn 



Mo 



+ ••• + 



dfj. 



MO 



X + 



where 



2n 



)fc k 

n k 
2n+l 



1,. 

n, 



1, 



(4.24) 



p{fi) = -LlYlifi-rUr) and p(^) = JJ (^ 



rrir 



r=2 



r=2 



in case of equations (HD) or ( ^W]) respectively. Here (pi, . . . ,(j)n are constant phases. This system 
contains n — 1 independent holomorphic differentials defined on the genus g = n—1 Riemann surface 
{w'^ = p{fJ^)}, which can be identified with the curve C described above. However, in contrast to the 
system ( ^.31) , in case of a polynomial p(/i) of odd order which corresponds to equation ( |HD| ), the 
last equation in (4.24) contains a meromorphic differential of the second kind having a double pole 
at the infinite point oo on C In case of a polynomial p{fi) of even order corresponding to equation 
( |SW| ) , the last equation includes a meromorphic differential of the third kind with a pair of simple 
poles at the infinite points oo_, oo+ on C. 

According to Clebsch and Gordon [1866] and Gavrilov [1999], in the odd order case, such a sys- 
tem describes a well defined and invertible mapping of the symmetric product C^9~^^^ to Jac(C, oo), 
the generalized Jacobian of the curve C with one distinguished point at oo. The set Jac(C,oo) is 
a noncompact algebraic variety which is topologically equivalent to the product Jac(C) x C. To 
describe this case we introduce a normalized differential of second kind having a double pole at oo. 



— ILqp^ dp 



+ ^dk 

k=i 



n 



1, 



(4.25) 



where LOk are the normalized holomorphic differentials specified in (|4.4|) , dk are normalizing con- 
stants such that all ^-periods of on C are zeros. Then the last equation in ( [4.24 ) implies 
that 



i=i •'w 



Z = V — ILqx + {d, Dt) + const, 

d = (di, . . . , dn-l)'^, t = {tn, . . . , ^2)"^, 



(4.26) 



where D is an (n—1) x (n — 1) normalizing matrix defined in ( |4.1lD . 

Since 00 now is a pole of , we choose the basepoint Pq = (/^O; ''^o) to be a finite Weierstrass 
point on C. For concreteness we choose Pq = (m2n,0). Applying the residue theorem to the 
generalized theta-function associated with Jac(C,oo) we solve the inversion problem ( [4.24 ) and 
find the following expression 



1=1 



2 2Zdv9[A + rj2n]{z) - dl9[A + r,2n]{z) 
Ol — Z H 



9[A + mn]iz) 
^Lqx + {d, Dt) + Zo 

Zq, zq = const, Ci = 



(4.27) 



z = Dt + zo£ C''-\ 
9 

pOk+ m2n, 



k=l 
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where the half- integer characteristic r/2n labels the point {m2n, 0), the vector V = {Dig, . . . , DggY^ 



is specified in (|4.4D , and the constant Ci contains the sum of integrals along the canonical cycles 



Ai, . . . , Ag on C. Notice that in the above formula dy = dt2- Expression ( 4.27] ) is meromorphic in 



X and ti, . . . , tn-i and can be regarded as a generalization of the Matveev-Its formula to the case 
of the noncompact variety Jac(C, oo). 

In the case of an even order curve C, corresponding to finite-gap peaked solutions of equation 
(|SW| ), system ( [1.24 ) defines a mapping of the symmetric product C^^^^^ to the generalized Jacobian 



Jac(C,oo-|-) which is topologically equivalent to the product Jac(C) x C*. As above, we set Pq to 
be the last Weierstrass point (m2n+i, 0) and introduce the normalized differential of the third kind 
having a pair of simple poles at oo_, oo+ G C, as well as the corresponding variable Z: 

where (di, . . . ,dg) = d are chosen such that all the ^-periods of Ooo± are zeros. Then, applying 
the residue theorem to the generalized theta-function associated with the Jac(C,oo-i-) yields 

5^ /i. + m = const ^ ^ (4.29) 

where, in view of ( [4.28| ), 



\J oo 



Z = x+{d,Dt) + Zo, z = Dt + zoeC^, 
Oi, . . . , / Og] G C^, Zq, zq = const. (4.30) 



Remark. According to the formula ( 2.22| ), expressions (4.27) and ( [4.29| ) describe formal solutions 



to equations ( |HD| ) and ( |SW| ) respectively. However, while treating these solutions, one needs to 
take into account the reflection phenomenon described in Theorem 4.1. Namely, when a certain 
variable //j passes zero, the point Pi = [ni, ypipi)) jumps from one sheet of the Riemann surface C 
to another or, in other words, from the pole of the differential of the third kind JIq to another 
pole Q_. Therefore, the above expressions do not provide global solutions to the equations. Instead, 
the following theorem holds. 



Theorem 4.3 1). The time- dependent finite-gap peaked solution U{x,t) of l \HL\ ) consists of an 
infinite number of pieces in M" = (ti, . . . ,tn-i,x) described by meromorphic functions 

TT ( ,\ n 72 ^ 2ZMdv9[A + mn]{zN) - ^^^[A + r,2n]{zN) , , ^ ^ 

Un[x, t) = Ci - H — — — , iV G ^, 

6^[A + ri2n\{ZN) 

zn = Dt + Nq + zq, Zj^ = \/^^Lqx + {d, zm) + Nh + Zq, Zq, zq = const, 

rQ+ , ^ / rQ+ fQ+ \ 

t = (ti,...,t„_i), h= ng\ q=[ wi,...,/ idg) , (4.31) 

JQ- \JQ- JQ- J 



where Ci is the constant specified in (J72^. 



For a fixed N the corresponding piece UN{x,t) is bounded by nonintersecting surfaces Sn-i and 
in M" given by equations 

Sn = {x=PN{t)}, PN{t)= J-^ {dvloge[A + r]2n]{zN+q/2)-{d,ZN)-Nh). (4.32) 
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The adjacent pieces Ui\i{x,t) and Ui\f+i{x,t) are thus glued to each other along Sn, where 



UipNit),t) = Cr-dl log0[A + r?2„](^jv + 9/2). (4.33) 

2). The finite-gap peaked solution U{x, t) of ( jSW^ consists of an infinite number of pieces given 
by meromorphic functions 

Un{x, t) = const ^ J', . ^ ^^ NeZ, (4.34) 



Z]\f = Dt + qN + zq, Zj\f = X + {d, zn) + Nh + Z| 



0, 



t = {tn, . . . ,t2), Zo, = const, h= $7 



''00± 1 



where the vector q is described in ^-3(1) . The piece UN{x,t) is bounded by peak surfaces Sn-i and 
Sn defined as follows 

Sn = {x =PN{t)}, pjv(t) = const -log ^j^j^^^ ^|^(oV (^•^^) 

t)[/\\[ZN + q + g/2j 

The adjacent pieces UN{x,t) and UN+i{x,t) are glued together along Sjy, where 

U{pN{t),t) = const - dv log ^i^ji^^'^j . (4.36) 

t)[/\\[zN + q) 

Notice that along the peak surfaces, the solutions described in 1) and 2) have discontinuous 
partial derivatives with respect to x and ti, . . . , t^-i- 

Remark. By fixing all the times but t^ in the above expressions, one obtains 2-dimensional 
piecewise solutions Uj\[ix,tf:), whereas the corresponding sections of Sj\[,<Sn C {x,t) = M"' describe 
peak lines in (x, tfc)-plane. As follows from ( 4.32| ) and ( 4.35| ), the motion of the A^-th peak PN{tk) 



along X-axis is described by a sum of a linear function in tj. and a quasi-periodic one. The latter 
function becomes periodic in the case g = 1. 

Finally, after fixing all the times without exception, expressions ( [1.31 ) and ( 1.34| ) provide pieces 



of the stationary finite-gap peaked solution already described in Theorem 4.2. 

Proof of Theorem 4.3. According to Theorem 4.2, the profiles of finite-gap peaked solutions are 
associated with geodesic ellipsoidal billiards and billiards in the field of a Hooke potential. An 
impact point on the boundary of a billiard trajectory corresponds to a peak of the profile U{x, to), 
and this happens when one of the passes zero. Hence, the solution ( |4.27 ) is valid until one 



of the points Pi, . . . ,P„ on C coincides with Q„ or 2+, the poles of the differential J^o in ( |4.5| 



Putting, for example, P„, = Q+ {^n = 0) in (4.24), one arrives at the following relations involving 



the normalized differentials defined in ( [4.4]) and ( 4.25| ), 

g 

' 'Po 



yZ ^k = Zk- qkl'i', k = l,...,g, (4.37) 

9 \ fQ+ / 9 

i=l " ^0 \ k=l / - ^0 \ 



YZ ( f^^^ - E d'^^'^) = ^Lqx - j^^ I 0« - dkuA , (4.38) 
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where Pq = (?n2n, 0). Notice that equations ( [4.37 ) form a closed system for the variables fii, . . . , fin-i 
and describe the standard Abel-Jacobi mapping C^^^ ^Jac(C). Hence, first symmetric polynomial 
has the following standard form in terms of theta-functions in the odd order case 



9 / 

Ml H h/in-i = ci - dv\ogd[/^ + 7]2n]{z -q/2), '^1 = XI r 



(4.39) 



On the other hand, the equation (4.38) implies that at a peak point the coordinate x becomes a 
function of z and therefore of t: x = poit). In the odd order case, this equation contains a sum of 
Abelian integrals of the second kind, so-called Abelian transcendent. By making use of the following 
standard expression for the normalized transcendent (Clebsch and Gordon [1866]) 



1=1 ■'M 

from (|3^) and ( pSD we find 
1 



Po 



Po{t) 



-IL, 







{dvloge[A + r,2n]{z-q/2)- {d,z)+h/2), 



Q+ 



(4.40) 



Using the trace formula for the solution U{po{t),t) = fii + - ■ ■ + and expression ([4.39| ) it follows 
that the equation x = pq (t) determines a surface Sq in C" along which the solution U has a peak. 
Now setting in ( |4.24D = Q- and taking into account (4^), ( 4.25| ) and 



Q- 



Po 



Po 



Po 



Q+ 



CO 



Po 



we obtain an expression for another peak surface Si determined by the equation {x = pi{t)} with 

1 



Piit) 



-ILf 



{dv log e [A + r]2n] {z + q/2)-{d,z)-h/2), 



(4.41) 



along which 

Uipi{t),t)=fii + ---+fin-i=Ci-dl.loge[A + r]2n]{z + q/2). (4.42) 

Under the reality condition, the surfaces Sq and Si do not intersect and therefore determine a 
connected domain in C" = {x,t) where the solution ( 4.27 ) is applicable. We denote this piece of 
solution as Ui{x, t). As follows from ( [4.40| ) and (4.41) Si is obtained from Sq by changing the phase 
as follows 



Z— + z ^ z + q that is x — > x + 



1 



-IL, 



-{h-{d,Dt)), t^t + D~\, 



(4.43) 



In addition, according to ( f4.42| ) and ( f4.39| ) at any two points on Sq and Si which are equivalent 
modulo the shift, Ui{x,t) takes the same values: 



Uiiqiit),t) = Ui{qoit) + 



l=.-ih-{d,Dt)),t + D-\) . 



(4.44) 



Now let us define function U2{x, t) = Ui{x + (h — (d, Dt)) / {\J —\Lq) ,t + D ^q), which is also a local 
solution to (HD). In view of i^^A^), Ui and U2 take the same values along Si, which ensures a 
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correct gluing of two pieces together. By using iteration with respect to both positive and negative 
N's, we construct complete sequence of peak surfaces and obtain formulae given in part 1) of the 
theorem. 

Similarly, solution ( [1.29 ) of ( ^W| ) is valid until one of the points Pi, . . . , Pn on C coincides with 



Q- or Q_|_, the poles of flo. Setting P„ = in (|4.24[ ) for the case of an even order curve C, and 



using (O) and (4.2g|) yields 



/ uJk = Zk - qkf^, k = l,...,n-l, (4.45) 
^ \ J^ooi - ^ dkCdk j = x- ( Jloo± - Yl j , (4.46) 



where Pq = (m2n+i)0). Inverting ([4.45| ) results in the following expression for a symmetric poly- 
nomial (see e.g., Clebsch and Gordon [1866]) 

m + • • • + ^^n-^ = const - dy log (4.47) 

After applying the theta-functional formula for the normalized transcendent of the third kind 
(Clebsch and Gordon [1866]), 

V- P 1 0[A](s-g) [^'^ 

gy^^^ f^oo,= const -log ^ = E^^' ^ = --1' 

from ( [4.46| ) and ( [1.45| ) we obtain 

X = poit) = const - log - d) + h/2. (4.48) 

By choosing P„ = in ( [4.24| ), one arrives at the expressions (l44^) and (^ ) with q/2,h/2 
replaced by —q/2, —h/2. Then, following similar arguments and applying induction, the piecewise 
solution of part 2) is constructed. 

We emphasize that although the different pieces Un{x, t) of the solution are obtained by iterative 
shifting the phases z, Z by the same vector, the pieces Uiy{x,to) of the solution (to being fixed) are 
all distinct because the shift occurs in both x- and t-directions. 

Remark. If we omit the reality condition above, then the hypersurfaces Sj\f, Sn in C" intersect, 
bounding a set of n-dimensional domains adjacent to each other in a rather complicated manner. 
Then the procedure of gluing different pieces of the functions UN{x,t) meromorphic inside each 
domain cannot be defined uniquely. As a result, the generic complex solution U{x,t) branches 
along the peak surfaces. 



5 Kinematics of peaks. 



Now we obtain expressions for the velocity of A^-th peak PAr(t) of the piecewise solution of (HD) 
with respect to time tk- As was shown above, the solution has a peak when one of the /i- variables 
passes zero implying that P„ = Q_ or P„ = Q i . 
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Theorem 5.1 Let yi, . . . , y-n-i denote the ^-coordinates of the points Pi, ... , Pn-i at the moment 
in time when one of the ^-variables passes zero. The following system of equations holds: 

dpN{t) ^ . s 

— = -Sfc-i(yi, . . . ,yn-i), (5.1) 

Otk 



where Syfc„i is the symmetric function defined in (2.25). In particular, we have 

= yi + • • • + yn-i = U{pN{t),t), (5.2) 

Ot2 

i.e., the t2-velocity of the peak coincides with its height. 

Proof. After applying limit mi — > 0, the equations (|2.12| ) and (|2.23| ) for the derivatives of take 
the form 



dx {Hn - /il) • ■ ■ (/U„ - fin-l) ' 



(5.3) 



-— = Sfc_i(/Xi,...,/X„_i ^ -. 5.4) 

otk (/i„ - /Zi) • • • (/i„ - /in_i) 

On the other hand, along the peak line {x = PN{tk)}, we have 

^Mp^(t.),t.) = + ^ = 0, 



which, in view of ( |5.3D and (|5.4D and after setting /i„ = 0, yields 

— -TT- + 5]fc-i(yi, . . . ,yn-i) = 0. 

1^1- fJ-n-l \ otk J 

Since p{0) 7^ and ni = yi, . . . , /U„_i = yn~i are finite, the latter relation gives ( |5.lD . 



Remark. The relations ( ^.1|) c an be also found by using direct differentiation of the expression 
for the A^-th peak surface ( |4.35| ) with respect to tk. Namely, putting without loss of generality 
N = 0, and taking into account dy = dt2, we write 

dpo{t) 



dtk 



dt.dt, iog0[A + mn] . (5.5) 



According to Mumford [1983], in case of odd order hyperelliptic curves, this gives a theta- functional 
expression for the coefficient in front of A""*^ in the polynomial (A — ^1) • • • (A — /.tg) which coincides 
with Sfc_i(yi,...,?/„_i). 



6 The dynamics of peaks and weak solutions. 

Expression ( |5.2D states that for equations in the hierarchies of (HD) or (SW), every peak in the 
solution profile move with velocity determined by the local value of the solution. In this section, we 
derive this property without recourse to tools related to the complete integrability of the evolution 
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equation. Thus, this property of peak motion can hold in general for equations that admit piecewise- 
smooth weak solutions, with jumps in the first spatial derivative at isolated points in the solution's 
support. In this case, the derivative discontinuity can be viewed as a 'shock' in the appropriate 
weak form of the evolution equation. 

We will take the weak form of the equation (HD) or (SW) to be 



/V^(.,*)^V(..*)...* = 0, 



(6.1) 



where the equality is satisfied for all test functions (j){x, t) is C°° with compact support in a domain 

in the (x, t) plane. Here = 
Y{x,t) = {Vi,V2) is defined by 



the dot denotes M inner product, and the vector function 



for equation (HD) and 
Vi 

V2 



Vi = 

V2 = 

= 5, 



y\ {Ul - 2kU) dy 



(6.2) 



2 4 



{2U^ + U^-2kU) dy 



(6.3) 



for equation (SW), respectively. We will look for jump conditions satisfied by the solutions of 
equation (|6.1|). If the jump discontinuities are isolated, by adjusting the support of the test functions 
(l){x, t) we only need to consider the case of a single discontinuity. 

Let us suppose that the function U{x,t) is infinitely differentiable almost everywhere in 0, 
except along the curve x = q{t) where the first derivative Ux has a discontinuity. If we partition 
the domain into O = U ^2 by cutting along the portion of the discontinuity curve x — q{t) = 
in 0, the divergence theorem and the choice of test functions (p{x, t) vanishing on the boundary 557 
allow us to write equation (|6.1| ) as 







dx dtVcj) -V 



dxdt<l)V-V+ dxdt(l)V-V+ / dl(l)n-[V 



(6.4) 



Here the unit vector n is directed along the normal [—q, 1] to the discontinuity curve dUi D 8^2 in 
Q, and [V]^ denotes the jump of the vector V across this curve. 



[V]! = lim V(x,t) 



lim V(x, t) 

x-*q{t)- 



By the arbitrariness of (^(x, t), each integrand term on the right-hand side of ( |6.4D has to vanish 
separately. Thus, from the first two terms. 



V-V = 0, 



or 



dVi dV2 
dt dx 



0. 



(6.5) 



in Oi or ^2, where U{x.t) is smooth. This smoothness and zero divergence condition, by the 
definition ( |6^ ) or ( |6.3| ) for (HD) or (SW) respectively, imply that ?7(x, t) is a solution of these 
equations in Qi or Q.2- For instance, (|6.5[) becomes 



Uxt + dx 



x-y\iU^- 2kU) dy 
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which is the integrated form of the Harry-Dym equation (HD). 

The last (jump) condition in (^), n • [V]l =0 along implies 

q[Vi]- = [y2]t. (6.6) 

The left-hand side of this expression is simply 

q[Vi]^ = Qpxt . (6.7) 



As to the right-hand side, the second (integral) term in the definitions ( |6.2D or ( |6.3D of V2{x,t) is a 
continuous function of x, as the integral wipes out the discontinuity sgn(x — y) as well as additional 
ones that Uy might have. Hence the integral terms do not contribute to the right hand side of 
(|6.6|) . The jump of V2{x,t) across the discontinuity curve x = q{t) then reduces to 

[V2]t = l[{U\)]^_ = U{q,t)[U.]t. (6.8) 

If 

[U,]+ = U,{q+,t)-U,{q~,t)^0 

equations (|6.7|) and (|6.8|) yield 



q=Uiq,t), (6.9) 

i.e., the location of the discontinuity (shock) in the Ux moves at the local speed U{q,t). We have 
then proved the following 



Theorem 6.1 Let U{x,t) be a solution of equation (6.1), with the vector V(a;,t) defined in terms 
ofU{x,t) by the nonlinear, nonlocal operators l \6. ^) and (6^) respectively for equations (HD) and 
(SW). Let U{x,t) be a smooth function of {x,t) in the domain Q C R^, except along the curve 
X = q{t), where U is continuous while the first derivative Ux has a jump discontinuity (peak) 
C/(g"'",t) 7^ U{q~,t). Then U{x,t) is a solution of equations (HD) and (SW) in each domain ill 
and O2 in which the curve x = q{t) partitions Q, and the location of the peak q{t) moves with 
velocity equal to its height, q = U{q,t). 



Conclusions. In this paper, profiles of the weak finite-gap piecewise-smooth solutions of the 
integrable nonlinear equations of shallow water and Dym type are linked to billiard dynamical 
systems and geodesic flows with reflections described in terms of finite dimensional Hamiltonian 
systems on Riemann surfaces. After reducing the solution of these systems to that of a nonstandard 
Jacobi inversion problem, solutions are found by introducing new parametrizations. The extension 
of the algebraic-geometric method for nonlinear integrable PDE's given in this paper leads to 
description of piecewise-smooth weak solutions of a class of A^-component systems of nonlinear 
evolution equations and its associated energy dependent Schrodinger operators. 
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